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I. INTRODUCTION 

The study of electronic transport of Dirac fermions in 
graphene is still one of the focuses in condensed mat- 
ter physics. Many theoretical works have been carried 
out since the experiments of graphene were realizedi 2 -^ 
Most of these studies (including the earlier ones in ad- 
vance of the experiment )!£i£iL&2Jii are based on the 
model of zero-range scatters or constant impurity scat- 
tering potentials in graphene. However, it has been 
found that the charged impurities with screened Coulomb 
potential o 11 ' 12 : 13 are responsible for the observed car- 
rier density dependence of the electric conductivity of 
graphene. 2 Therefore, a satisfactory transport theory for 
Dirac ferions in graphene should take into account the 
charged impurity-electron scatterings. 

Recently, we have developed a formalism for the elec- 
tric conductivity of Dirac fermions under the charged im- 
purity scatterings within the self-consistent Born approx- 
imation (SCBA).— The formalism has been extended to 
study the Hall conductivity of graphene.— Our results 
show that the experimentally measured electric conduc- 
tivity and the inverse Hall coefficient can both be suc- 
cessfully explained in terms of the carrier scatterings off 
charged impurities. Based on the same scattering mech- 
anism, we have studied the weak localization (WL) prob- 
lem in graphene^ by considering only the conventional 
maximum crossing diagram.— In that work, the reason 
why the WL in graphene has not been detected in the 
experiments was given. In the present paper, we go be- 
yond our previous approach^ by including two additional 
diagrams as suggested in Ref. la, and reexamine the WL 
effect in graphene taking into account realistic situations 
of experiments. The major part of our effort is to solve 
the Cooperon propagator under the finite-range scatters, 
which is quite different from the one under the zero-range 
scatterings as studied in the existing works 
The detailed formalism for the Cooperon propagator is 
presented. With the Cooperon propagator obtained in 
a weak magnetic field, we calculate the magnetoconduc- 
tivity of graphene and compare the results with recent 



experimental measurements i 16 ' 17 

The electron system in graphene at low energy excita- 
tions can be viewed as massless Dirac fermions ^ 18 ' 19 ! 20 
as being confirmed by recent experiments.—^ Using the 
Pauli matrices ex's and r's to coordinate the electrons in 
the two sublattices (a and b) of the honeycomb lattice 
and two valleys (1 and 2) in the first Brillouin zone, re- 
spectively, and suppressing the spin indices for briefness, 
the Hamiltonian of the system is given by 
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where ^ = {c\ al , c\ b} , c\ b2 , c\ a2 ) is the fermion opera- 
tor, the momentum k is measured from the center of each 
valley, v (~ 5.86 eVA) is the velocity of electrons, V is 
the volume of system, and Vi(q) is the finite-range im- 
purity potential. From our previous result* 2 ^ the cutoff 
for k summation is about k c w ir/3a (with a the lattice 
constant), within which the electrons can be regarded as 
Dirac particles. For charged scatters, Vi(q) is given by, 



Vi(q) = 
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where rii(—q) is the Fourier component of the impurity 
density, vo(q) and v\ are respectively the intravalley and 
intervalley impurity scattering potentials, and Q is a vec- 
tor from the center of valley 2 to that of the valley 1 
[Fig. IDja)]. In Appendix, we detail the discussion on 
this impurity potential. Here, all the momenta are un- 
derstood as vectors. We will adopt the charged impu- 
rity potential as the Thomas- Fermi (TF) type, vo(q) = 
2ire 2 exp(— qRi)/(q + <?tf), where Ri is the distance of 
the impurity from the graphene layer, qTF = ^kpe 2 /ve 
is the TF wavenumber, kp is Fermi wavenumber, and 
e ~ 3 is the effective dielectric constant. (In terms of the 
carrier density n, kp is given by kp — y/im. We will use 
the parameter S defined as the doped carriers per carbon 
atom to denote the carrier concentration. It is related 
to n via 5 = \/3a 2 n/4.) The intervalley scattering po- 
tential is defined as v\ — vq(Q) with Q = 47r/3a as the 
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where uj = Eo(fc, w) with fi the chemical potential, 

hf. = vk + S c (fc, u>), k is the unit vector in k direction, 
and the frequency uj is understood as a complex quantity 
with infinitesimal small imaginary part. 

According to the SCBA, the current vertex is made of 
the ladder diagrams as shown in Fig. Hljc). The current 
vertex vT x {k, u>i, UJ2) is expanded as 



r x (fc,wi,w 2 ) = ^2y j (k,uj 1 ,uj 2 )A*(k) 



(5) 



FIG. 1: (color online) (a) Brillouin zone and the two Dirac- 
cone valleys, (b) Self-consistent Born approximation for the 
self-energy, (c) Current vertex with impurity insertions, (d) 
Current-current correlation function. The solid line with ar- 
row is the Green function. The dashed line is the effective 
impurity potential. 



distance between the two nearest Dirac points. In order 
to fit the experimental data, the impurity concentration 
is chosen as ni = 1.15 x 10~ 3 a~ 2 M- The distances Ri 
may be within the range of a few lattice constants. Be- 
cause the typical momentum transfer q due to intravalley 
scattering is of the order of fcp, qRi << 1 for the doping 
levels considered here, the electric conductivity a does 
not sensitively depend on Ri. Comparing to vo(q), v\ 
is very small unless some of the impurities are close to 
the graphene layer. For the case of the impurities on the 
layer, we have v\ = 2ite 2 /(Q + (/tf)- In our following 
discussion, we set Ri = a and v\ = a2ire 2 /(Q + <7tf) 
with a an adjustable constant. 

To theoretically study the electronic transport in 
graphene, we need to start with the Green function 
theory i 22 i 23 ' 24 i 25 ' 26 We here briefly review the main re- 
sults for the SCBA and the current vertex correc- 
tion. Under the SCBA [Fig. ffl^b)], the Green function 
G(k,Lu) = [uj + fx — vk ■ <7T Z — E(fc,w)] _1 and the self- 
energy uj) = So(fc, uj) + S c (fc, u>)k ■ Bt z of the single 
particles are determined by coupled integral equations: 13 



where A^k) = t z <j x , Af(k) = a x B ■ k, A 2 (k) = B ■ k<r x , 
A^(k) = t z B ■ ka x S ■ k, and yj(k,u>i,U!2) are determined 
by four-coupled integral equations?^ The ^-direction 
current-current correlation function [Fig. [U^d)] is ob- 
tained as 



P(u>i,w 2 ) = 
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gTr[Af (fc)G(fc, LOx)A^(k)G(k, uj 2 ) 



and 4> is the angle of the vector fc, for w's (uj\ and ui 2 ) 
= uj ± iO = cj ± . The detailed derivations of T x (k, uji,uj 2 ) 
and P(uji,lu 2 ) can be found in Ref. IT3L and will not be 
repeated here. 



II. COOPERON 

A major task of theoretically studying the quantum 
interference effect in the electronic transport is to find out 
the solution to the Cooperon propagator in the presence 
of charged impurities. Although this problem is discussed 
in Ref. [TEI, the detailed derivation has not been given 
there. In subsection A, we solve the Cooperon propagator 
of Dirac fermions in the absence of external magnetic field 
B. In subsection B, we present the solution to the case 
with the existence of the magnetic field. 
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E (fc,w) = ^^(Ifc-fe'D+^SoC^) (3) 

k' 

E c (fc,w) = ^^ 2 (|fc-fc'|) ffc (fc>)fc-fc' (4) 

k' 



The Cooperon propagator C^^ 4 Q4 (fc, k 1 , q, uj) obeys 
the Bethe-Salpeter 16 x 16 matrix equation represented in 
Fig. 2(a). Here, the superscripts j's denote the valley in- 
dices, and the subscripts a's correspond to the sublattice 
indices. To explicitly write out the equation of Fig. 2(a), 
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we here give the simpler one for C^^^ ai (k, k', q, w) 
that starts from the single impurity line [the dashed 
line with a cross in Fig. 2(a)], using the convention 
^mct 2 = ^hj2^a 1 a 2 an d 3 {&) as the conjugate valley (site) 
of j (a): 



(b) 
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FIG. 2: (color online) (a) Cooperon propagator in momentum 
space and (b) in real space. 
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The first term in the first line in right hand side of 
Eq. {B{ is due to the intravalley scatterings, while 
the second term comes from the intervalley scatterings. 
Sj ± j 2 means that when a particle is scattered to val- 
ley ji, another particle should be scattered to valley j\ 
because of the constraint for the intervalley scattering 
(ni(Q)rii(Q')vf) /V = rviv\§Qi-Q. It also ensures that 
the total momentum (vanishing small under considera- 
tion) of the Cooperon is unchanged. The second and 
third lines are the processes of Cooperon propagating 
after the intravalley and intervalley scatterings, respec- 
tively. The equation for the Cooperon propagator C is 
obtained by subtracting the single impurity line from C 
in Eq. ©. 

The form of Eq. ^ seems rather miscellaneous. It 
may be simplified by classifying it with good quantum 
number of the Cooperon. To do this, we note that the 
elements of the coefficient matrix of C in Eq. ^ are 
arranged according to the indices (superscripts and sub- 
scripts) of the Green functions. Since the Green function 
G(A;,Cl>) are composed by the unit matrix and t z <j ■ k, we 
then look for all the operators that commute with t z g ■ k. 
Recently, McCann el al& have introduced the operators 
of isospin E's and pseudospin A's, 

E() = T C7o, S X =T 2 (T X , E y =T z (T y , £ Z =T (7 Z , 

Ao =T a , A x ~T x cr z , K y —T y a z , A z — r z a . 

Clearly, A's commute with S's and t z g ■ k, and are 
conserving operations for the Cooperon. Therefore, we 



transform the Cooperon from the valley-sublattice space 
into the isospin-pseudospin space according to McCann 
et ai, 

C U ' = - \ (M l V 102 niifaiaji CAtl'\\3&ja (<7\ 

SS' ^ / , l lvJ sMiQ2 (11020314^ S r )<X4(X3 \' I 

where M\ = E^E^A^A/. We will hereafter occasionally 
use the indices 0,1,2,3 or 0,x,y,z to label I and s. In the 
isospin-pseudospin space, the single impurity line is given 
by W l £ = W l s S l £, with 

W l s (\k - fci|) = mvl(\k - fell) + n lV l{8 w - S lz )(-l) s , 

which is the transform of the first line in the right hand 
side of Eq. ([6]). The result of second+third lines in right 
hand side of Eq. is transformed to 

±^[W l (\k-h\)h(ki,q)} SSl Cl l L'(ki,k',q) 

fci,si 

where h(k\,q) is a matrix defined in the isospin space 
with the element given by 

/W(fci,<?) = Tr[G(-fc+ w+)E,G(-A;r,w-)Et,]/4 (8) 

and kf = ki ± q/2. With these results, we obtain the 
equation for C l J s , , 

C l Uk,k',q) = y ^nL 1 (fc,fc 1 ,g)«(|fc 1 -fe'|)C' 

fel,Sl 

+ C»' s ,(fc 1 ,fe',< Z )] (9) 



4 



where tl l (k, ki,q) — W l (\k — ki\)h{ki,q). Here, the ar- 
gument oj of C l j s , and II 1 has been suppressed for brief- 
ness. From Eq. (|9|). it is seen that the pseudospin of the 
Cooperon is indeed conserved during the impurity scat- 
terings, C l J s , — C l ss ,5w. We then need to deal with C l ss ,. 
Thus, the original 16 x 16 matrix equation is separated 
into four 4x4 ones, each of them corresponding to a def- 
inite pseudospin I. In the isospin space, the Cooperon of 
a pseudospin I is a 4 x 4 matrix denoted as C l . 

To solve Eq. ([9]), we use the standard method that 
expands C l in terms of the eigenfunctions W n (k,q) of 
n'(fc s fci,g): 

C\k, k', q)=J2 c l n (q)<(k, q)^(k', q), (10) 

n 

where c l n (q) are constants and 

^J2 Rl ( k ' k ^l) iS[ n(ki,q)=X l n (q)^ n (k,q) (11) 
fei 

with X l n (q) as the eigenvalue. Here, ^ l n (k,q) is a column 
vector with four components in the isospin space since 
H l (k, ki, q) is a 4 x 4 matrix in this space. The constants 
c l n (q) are determined by substituting Eq. (fTQ|) into Eq. 
©. It is then seen that c l n (q) oc [1 - X^q)}' 1 . Therefore, 
the predominant contribution to C l comes from the state 
with the lowest |1 — A„(g)| ee |1 — X l (q)\ that can be 
vanishing small. We will here take into account only the 
state of the lowest |1 — A'(g)| for each I. 

Firstly, we consider the case of I = and q = 0. A 
solution can be explicitly obtained as A°(0) = 1, and 
*°(fc,0) EE with 

**(fc) = [A o (fc,w),-A c (fc,w)cos0,-A c (fc,w)sin</>,O] 

where ^*(fc) is the transpose of ^(fc), Ao(k,oj) — 
ImEo(fc,w _ ), and A c (fc,cj) = ImS c (fc.w"). The four 
components of ^(fc) correspond to s — 0,x,y,z respec- 
tively. The solution of A°(0) = 1 is the most important 
one which gives rise to the diverging contribution to the 
Cooperon. One may check this result with the help of 
Eqs. ([3|) and (T4|. Actually, the above solution is just a 
consequence of the Ward identity (under the SCBA for 
the self-energy): 

Tx6S(k,u-) = ^^{V^k - k')1mG{k' ,u>-)Vi{k' - k)) 

k' 

where (• • • ) means the average over the impurity dis- 
tributions [Fig. HKb)]. There are three non-vanishing 
components in 'J'(fc) because of the finite-range impurity 
scatterings. For the zero-range potential, only the first 
component of survives and is a constant. One then 
needs to solve a scalar equation instead of the matrix 
integral equation. 

For finite but small q, by expanding fl (A;, fc', g) to 
second order in q and regarding the difference from 
II (fc, fc',0) as a small deviation, we then solve the 



problem by perturbation method. Since expanding 
IL°(k,k' ,q) [equivalent to expanding h(k',q)] is an ele- 
mentary manipulation but tedious [because 9 elements 
in h(k',q) need to be expanded], we here just present 
the result. For / ^ 0, the difference between Ir and II 
comes from the intervalley scattering term in W l . Sim- 
ilarly, we can treat this difference by perturbation. For 
all the cases, to the first order in the perturbation, we 
have 

« A'(0)-d ; g 2 , for q^O (12) 

with 

X l (0) = 1-(1 + 5 lz -S w )SX 

SX = n^(A )(A (M 2 + \9c\ 2 ) + 2 A c Re(g+g- ))/(*]*) 

(13) 

where (^l*) = Y.k^ (k)^{k)/V , (■■■) here means an 
average over the Fermi area, i.e., (A ) ee A (fc,o;)/V, 
9oc = 3o,c(fc,t^ ± ), and the arguments k and lo have been 
suppressed for clarity. The coefficient di is given by 

d t = (2D L A \g c /k\ 2 + (D l A + D c A c ){\g' \ 2 + \9' c \ 2 
-Re[ ff +(. 9o -" + go'/k) + g+(g-" + g-'/k)]} 
+ (D l A c + D c A )Re[2g+'g-> - g+ig^" + <fo 7*) 
-5o + (5 c ~" + 9c' /k- ff -/fc 2 )])/8(*|*> 

where D l = E k ' W o(\ k ~ k'\)A (k')/V , D c = 
Efe' Wl(\k-k'\)A c (k')cos9/V with the angle between 
k and k', and g' = dg(k,ui)/dk. To the 0th order, the 
eigenfunction ^S(k) is unchanged. We here consider only 
the case of small q since that is where QIC is significant. 

In Fig. 3, the eigenvalues X l (q) are shown as functions 
of q at electron doping concentration 8 = 2.0 x 10~ 4 
and the intervalley scattering parameter a = 0.7. The 
eigenvalue for I = x, y is degenerated. In the limit q — ► 0, 
only 1 — X°(q) approaches zero. As shown by Eq. (|13[) . 
the finite values 1 — A x (0) = SX (shown in the inset of 
Fig. 3) , and 1 — A 3 (0) = 2<5A sensitively depend on the 
intervalley scattering strength v\. At v\ = 0, A'(0) = 1 
is valid for alH. 1 — A'(0) is larger for stronger v±. On 
the other hand, the coefficients di are almost the same 
for all the pseudospins. 

The state for each I so obtained is of the lowest 1 1 — 
A'(<7)|. For the lowest state, c l n (q) is given by c l (q) — 
ci/[l-X l (q)} with 

01 = » E * t (*)n'(fc l fci,o)w"(|A 1 - fc'|)*(fc') 

where the (/-dependence of c; has been neglected because 
of the drastic behavior of 1/[1 — A (q)] at small q. The 
Cooperon is finally approximated as 

C l {k,k',q) = cc&{k)*\k>)l\l- X\q)]. (14) 
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FIG. 3: (color online) Eigenvalue X l (q) as function of q at 
electron doping concentration 8 — 2.0 x 10 -4 . The inset shows 
<5A as function of doping concentration 5. 

For the zero-range scatters, only the isospin-singlet 
C l m (k, k', q) survives and is independent of k and k 1 . 
In this case with v\ = and w = 0, by the one- 
band approximation, one obtains A (q) = 1 — Dq 2 /4Ao 
where D = v 2 /2A is the diffusion constant and A is 
the first component of 'F. The resultant Cooperon is 
AAoriiV 2 /Dq 2 consistent with that of Ref. d to the order 
of q~ 2 in q — > 0. 

B. B^0 

The Cooperon in the existence of weak magnetic field 
can be treated in a way parallel to Ref. [25|. Here, we 



outline the main steps. When a weak external magnetic 
field is applied to the system, we need to start with the 
description in real space because of its inhomogeneity. 
The kinetic part of the Hamiltonian Hq is 



H = / dfffi (r)vr z a ■ (p + A)ip(r) (15) 



where p is the momentum operator, and A is the vector 
potential. Here, we have used the units of e = c = 1 
(with — e as an electron charge and c the light velocity). 
With gauge transform 



ip(r) = tp(r) exp(— i / dl-A), 

J r 

A is eliminated from Hq. On viewing this gauge trans- 
form, for very weak magnetic field B, the Green function 
G(r,r',T — T r ) = —{Ttj}(r,T)ip*(i J ) T')) can be approxi- 
mated as2£ 



G{r,r',uj)^G(r-r',uj)exp{i dl-A), (16) 

where G(r — r , ui) is the Green function in the absence 
of the magnetic field. Here, the position r is understood 
as vector. 

In real space, the Bethe-Salpeter equation for the 
Cooperon is diagrammatically shown in Fig. 2(b). One 
can then write out it explicitly and do the same transform 
as from Eq. ((6]) to Eq. ([9]). The final matrix equation 
(in the isospin space) is given by 



C l {r 1 ,r 2 ,r 3 ,r 4 ) = W l (ri - r 2 )[h(r 1 , r 2 , r 3 , r 4 )W l (r 3 - r 4 ) + \ dx x / dx 2 h(r 1} r 2 , Xi, x 2 )C l (xi, x 2 , r 3 , r 4 )] (17) 



where W (r) is the real space representation of W (q), and the element of h(r\, r 2 , r 3 , r 4 ) is given by 

h S s'{ri,r 2 ,r 3 ,r 4 ) = ^Tr[G(n, r 3 , w+)E s G(r 2 , r 4 , a; - )!!*,]. 
Using the approximation given by Eq. (JTHJ) , we have 

h ss <{ri,r 2 ,r 3 ,r 4 ) rj h° ss , (n - r 3 , r 2 - r 4 ) exp(i / dl ■ A + i / d/ • A) 

« h° ss ,{R- R' + r —^- 1 R- R' - r — — )exp(i2 / df- A) (18) 
2 2 J fl 

where i? = (n + r 2 )/2, r = r% — r 2 , R' = (r 3 + r 4 )/2, r' = r 3 — r 4 , and h° is defined in the absence of A. The Fourier 
transform of h°(R - R' + r -=f, R-R' - ^f-) is given by 

h°(R -R'+ r —^,R — R' — r —^) = -^Y h(k, q ) e ^(R-R')+ik-(f-f) (ig) 
2 2 V ~* 

kq 
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where h(k,q) is defined by Eq. ([8]). The eigenvalue problem of Eq. (|17[) reads 



W l (n-r 2 ) / dr[ / dr' 2 h{r 1 ,r 2 y i y 2 )^{r' i y 2 )=E l ^{r 1 ,r 2 ). (20) 



where if)(ri,r 2 ) is a four-component vector in the isospin space, and E l is the eigenvalue. Using the coordinates R 
and r, we separate ip(ri,r 2 ) as ip(ri,r 2 ) — ^(R)^(r) where $(i?) is a scalar representing the motion of center of 
mass, and ^S(r) is a four-component vector meaning the relative motion of the Cooperon. Since the magnetic field is 
weak, only the large-scale motion of the center of mass is significantly affected; the magnetic field influence on the 
relative motion is negligible. Then, 'i'(r) can be considered as the real space representation of ^(k) given in Sec. II. 
By integrating out the relative motion, we get 

TfJ2 f dR'\ l (q)exp[iq- (R- R')+i2 [ dT- A]<S>(R') = E l $(R) (21) 

V q J J R 



with \\q) defined by Eq. (12]). Using \\q) » \ l (0)-d iq 2 
for small q, <3>(i?') = exp[(_R' - R) ■ V]$(i?), and carrying 
out the q-summation and -R'-integral, we obtain 

[A'(0) -dt(P + 2A) 2 ]<5>(R) = E l <f>(R), (22) 

where P = — iV is the momentum operator (of the center 
of mass) of the Cooperon. Using the Landau gauge A = 
(0, Bx 1 0), one has 

E l n = X l (0)-4d l B(n + 1/2), (23) 

and $„(i?) being the wavefunction of the corresponding 
Landau state. The degeneracy of each level is g — BV/ir. 
In real space, the Cooperon is expressed as 

&(n, r 2 , r 3 , r 4 ) = ff £ Q *(r)*V). 



III. QUANTUM INTERFERENCE 
CORRECTION AND 
MAGNETOCONDUCTIVITY 

The magnetoconductivity measured in 
experiment o 16 ^ 17 can be shown to stem from the 
quantum interference correction (QIC) to the electric 
conductivity. We here firstly discuss the QIC in the 
absence of the magnetic field. For the conventional 
electron systems, the main contribution to the QIC is 
theoretically given by the maximum crossing diagrams 
as shown by Fig. 4(a) i 23 ' 27 i 28 The works based on this 
diagram have also been done to understand the WL 
problem in graphene^i^ However, as it has been pointed 
out£ that there exist additional diagrams such as those 
in Fig. 4(b) which have the same order as Fig. 4(a). 
When the Fermi energy is not much larger than the 
level broadening due to the impurity scatterings, the 
contribution from other higher order diagrams may not 
be negligible^ In Fig. 5, we dipict the ratio of the 
Fermi-level broadening 7 and the Fermi energy Ep as 




-k+q/2 a)- k'+q/2 -k'+q/2 k+q/2 



FIG. 4: (color online) (a) Maximum-crossing-diagrams correc- 
tion to the conductivity, (b) Additional quantum interference 
correction to the conductivity. 



function of the carrier concentration 5 obtained for the 
impurity density m = 1.15 x 10~ 3 a~ 2 and the intervalley 
scattering parameter a — 0.7. At 8 < 1.0 x 10~ 4 , this 
ratio becomes larger than unity, the present approach 
based on the perturbation expansion my not be appli- 
cable. Our calculation for comparing with experimental 
measurements will be performed at S > 2 x 10~ 4 where 
•y/Ep < 0.4. It is important to notice that except for the 
ladder diagrams given by Fig. [ijd) which are generated 
from the SCBA, any theoretical schemes including the 
diagrams such as in Figs. 4(a) and 4(b) are not the 
conserving approximation. The QIC is meaningful and 
qualitative only when it gives small correction to the 
conductivity. In this sense, the results presented below 
are qualitative and applicable only when the carrier 
density is not too small. 

With the Cooperon C , we firstly calculate the QIC 
from Fig. 4(a) to the current-current correlation function 
5Pi(u>~ , cj + ). Because the vertex, the Green functions 
and the Cooperon are matrices, one needs to write out 
5P\(lo~ ,lu + ) carefully. For doing it, we start to work in 
the valley-sublattice space. According to the Feynman 
rule, we have 
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FIG. 5: (color online) Ratio between Fermi- level broadening 7 
and the Fermi energy Ef as function of carrier concentration 
8 calculated at impurity density m = 1.15 x 10 -3 a~ 2 and 
intervalley scattering parameter a = 0.7. 



where T4(fc,^i,w 2 ) = ui)T x (k, Wi, W2)G(k, 0J2) is 
the vertex connected with two Green functions. With 
the inverse transform of Eq. ([7]) 



we get 



r-ihhhii — _ \ " (jufl\\hh n 11 ' ( M 1 ' V' 

aiQ2tl3CK4 1 / j \ 8 J Ct20tl SS* \ s' J a 



334. 
3 Q 4 ' 



we obtain 



= ^^Tr[Z'(fc,c)C i (-fc ! fc, g )] (25) 



kql 



where Z l (k,uj) is a matrix with elements Z l ss , defined as 

Z\ s , = TT[Vt(k,U + ,LO~)M L s V x (-k,LU~,LO+)Ml?}. 

There is a simple relation, Z l (k, ui) = —Z Q (k, u) for / 7^ 0, 
because 

M l s = M°A h A;M S %* = -M°, for l = x,y,z 

and the operator A; commutes with G and T x . This re- 
sult means that the QIC by the pseudospin singlet (I = 0) 
is negative, while it is positive by the pseudospin triplets 
(I = x, y, z). Substituting the results given by Eqs. (fT5|) 
and ([T4"l) into Eq. (|25p and carrying out the g-integral, 



8P 1 (u>-,u,+) = fy2N l ^lnj 
^ di 1 



d t l-\ l (Q)+d iq 2 



, (26) 



f = -^y y £^WZ°(k,u)^-k), ( 27 ) 



where No = —1, N =x _ y _ z = 1, qo and gi are the lower 
and upper cutoffs of the g-integral. 

The lower cutoff qo is given by go = max (L~ n 1 , L~ x ) 
where Li„ is the length the electrons diffuse within an 
inelastic collision time Ti n and L the length scale of the 
system.^ At low doping 8 < 1.0 x 10~ 3 , r in due to 
the inter-electronic Coulomb interaction is estimated as 
Ti n ~ 0.462w/aT 2 (where T is the temperature) from 
the recent study of the interacting electrons in graphcne 
using renormalized-ring-diagram approximation. — L in is 
then given by L;„ = (v 2 TTi n /2) 1 / 2 where the elastic colli- 
sion time r is determined by the non QIC-corrected con- 
ductivity Co, t = hirao/vkpe 2 (with kp as the Fermi 
wavenumber).— For low carrier density, we find that Li n 
is about a few microns for 4 K < T < 20 K. On the other 
hand, the upper limit is qi = Lq 1 with Lq = vt as the 
length of mean free path. 

Similarly, we can obtain the contribution from dia- 
grams in Fig. HJb). Since the contribution from each 
of the two diagrams can be shown to be complex conju- 
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gate of each other, we then find out the formula from the 
first diagram, take the real part and multiply it by factor 



J 



2. The result is 



6P 2 (w-,w+) 



v 

V 3 



ReTr{C l ss ,{k,k' ,q)V,t(k,Lj + ,u)-)M l s ,G(-k,u)-)[V x {-k' ,uj + )M l s *G t {k' ,u)-)mvl{k - k') 



kk' qlss' 

+TiaiV x (-k\u}~ ,Lu + )M l s *G t (k' ,uj~)Tiainivl]}. 

I 



(28) 
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FIG. 6: (color online) Quantum interference correction to the 
conductivity as function of electron doping concentration 5 
for the intervalley scattering parameters a = 0.7 (solid line) 
and 0.5 (dashed line). 



Since the maximum crossing diagrams in Fig. SUb) in- 
clude the single impurity line, we here have 



C l ss , (k, k', q) = C\ s , (k, k', q) + W l s (k~ k')S s 



(29) 



Corresponding to the approximation for C l ss , (fc, k! , q), we 
also expand the matrix W l (k — k 1 ) in terms of the eigen- 
functions of the peudospin-singlet Cooperon and take 
only the ground-state contribution. Then, the (/-integral 
of C l ss ,(k, k', q) gives rise to a logarithm term and an ad- 
ditional term due to respectively C l and W l in the right 
hand side of Eq, $Z§§. With comparing to 8P 1 (0~,0 + ), 
<5P 2 (0\ 0+) varies in the range of -0.3 to -0.7 <5Pi(0-,0+) 
for carrier concentration S < 2 x 10~ 3 at low temperature. 
This feature is some what different from that of Ref. 6 
in which the ratio is -0.5 independent of 5 for zero-range 
impurities. 

With the current-current correlation function, one then 
calculates the conductivity a according to the Kubo 
formalism^ At very low temperatures, the correction to 
the conductivity is calculated by 5a — 5P(0~, + )/2n. 



The QIC to the electric conductivity as function of elec- 
tron doping concentration is shown in Fig. 6. For the 
intervalley scattering parameter a = 0.5, the QIC is pos- 
itive at low doping for the sample size L = lfim, im- 
plying the system may be delocalized. While at large 
doping the correction becomes negative. Physically, at 
small electron doping, the Fermi circle and the typical 
momentum transfer q (~ 2kp) are small and also the 
screening is weak, leading to stronger vo(q) than v±. For 
very weak intervalley scattering, all A'(0)'s (I ^ 0) close 
to 1 (see the gap SA as function of 5 in the inset of Fig. 
3), the QIC from the pseudospin-triplet channels has al- 
most the same magnitude. After one of I ^ is canceled 
by the I — channel, the net QIC is positive. This is 
consistent with the fact that Dirac fermions cannot be 
scattered to exactly the backwards direction in case of 
vi = and the WL is absent. On the other hand, with 
increasing electron doping, the strength of the intervalley 
scatterings becomes stronger, leading to the appearance 
of WL in large size samples. The strength of QIC can be 
weakened/quenched when the sample size < L;„ (about 
a few microns at low temperatures) as studied in the 
experiment. 2 ' 30 In large size samples, the WL effect is 
observable at finite carrier doping. The present result is 
qualitatively the same as of our previous studyi^ where 
the numerical value, however, is not accurate due to a 
numerical error. The present model of the impurity scat- 
terings and the results should be more realistic than the 
previous ones. 

We here concisely explain why pseudospin singlet 
Cooperon give rise to WL (negative QIC) but the pseu- 
dospin triplets result in anti-WL (positive QIC). As seen 
from the matrix M l s defined below Eq. ([7]), the pseu- 
dospins are actually associated with the matrices A^A;, 



A y A = 

AyA X = 
AyAy = 

AyA z = 



t 2 <7 3 , singlet 

— 1x3, triplet with I = x 

1, triplet with I = y 

iTi(73, triplet with I = z. 



In terms of the c-operator, ip t q / 2 -k^y^-i' l Pq/2+k annihilates 
a Cooperon of total momentum q and relative momen- 
tum fc with pseudospin I. It is seen that in the pseudospin 
singlet state the two particles are in different valleys and 
the parity is odd under the exchange of the valley indices. 



9 




FIG. 7: (color online) Real space representation of the quan- 
tum interference correction to the conductivity corresponding 
to Fig 2(a). 

The magnitude of the wave function of the pseudospin 
singlet Cooperon is large when the two particles occupy 
respectively the opposite momentum [defined respect to 
the origin of the Brillouin zone, see Fig. Ula)] states of 
the single particles. It implies a strong backward scat- 
tering for the electrons and thereby leads to WL. On the 
other hand, the pseudospin triplets are even under the 
valley exchange. For I = z, though the two particles are 
in different valleys, the wave function of the Cooperon 
is small when the two particles occupy respectively the 
opposite momentum states. In this case, the backward 
scattering is weakened, resulting in the increase of the 
conductivity. For / = x and y, the two particles occupy 
the states in the same valley and their total momentum is 
finite. The case corresponds to the final state of the scat- 



tered electrons is not in the backward direction, giving 
rise to a positive contribution to the conductivity. The 
anti-WL can be reduced when the backward scattering is 
strengthened by the intervalley scatterings. 

The QIC to the conductivity in the presence of mag- 
netic field is calculated according to Figs. 4(a) and 4(b). 
These diagrams can be redrawn in real space. For ex- 
ample, shown in Fig. 7 is the real space representation 
of Fig 4(a). To explicitly write out the expression, we 
note that the magnetic field effect on the vertex is negli- 
gible small. The strong dependence of B in QIC comes 
from the Cooperon due to the small denominator I — E n . 
Therefore, the vertex associated with x' in Fig. 7 con- 
nected with two Green functions after integrating over x' 
is given by 

V* (r 3 - ra, w+ w") = y X) V * 0J+ ' w ~) e ^ ' ^ ' 

fei 

where V x (ki, uj + , oj~) has appeared above Eq. (|25|) . Sim- 
ilarly, for the vertex associated with x in Fig. 7, one gets 

V x (r 4 -n,uj-,uj + ) = ^^2v x (k 2 ,u>-,u> + )e i ^-^. 

The QIC to the current-current correlation function from 
Fig. 7 is 



SP 1 (lj-,u; + ) = ^Y1 J dri '"J ^4Tr[Z'(r3-r 2 ,r4-r 1 , W )C'(r 1 ,r 2 ,r3,r 4 )] (30) 

with Z l ss ,(r 3 - r 2 ,r 4 ~ n,u) = Tr[F x *,(r 3 - r 2 , uj+, Lo~)M l s V x {r A - n, u~ , u + )M l s f]. Substituting the results for V x >, 
V x , and C l (ri, r 2 , r^, n) into Eq. ([30]) . setting ki, 2 — ±k + q/2 and neglecting the small g-dependence in V x (±k + 
q/2, ui, UJ2), then using the coordinates R, i?', r, and r', and integrating out the relative motions, we get 

Inq n ' In n 



By comparing this result with Eq. (j25|) . we see that the 
g-integral in Eq. (|25p is replaced with the summation 
over the Landau levels. By using the same g-cutoffs as 
in obtaining Eq. ([^rj]) . only those states of energy levels 
e l n = B(2n + I) (with E l n = X l (0) - 2dit l n ) in the range 
(<$/2,gf/2) (in units of a = v = 1) need to be summed 
up. Similarly, we can obtain the QIC 8P2 from Fig. 0Jb) 
or Eq. (|28p in the presence of the magnetic field. 

The magnetoconductivity is defined as Aa(B) = 
cr(B) — cr(0) where a is the corrected conductivity in- 
cluding the non-corrected one and the correction da. 

Shown in Fig. 8 are the present results for the magne- 
toconductivity Aa(B) and comparison with experiments. 
For comparing with the experiments, we note that the 



overall magnitude of Aa(B) varies largely from sample 
to sample (of the same carrier density) and from experi- 
ment to experiment ] 16 ' 17 The reason is that the interval- 
ley scatterings varies in samples and in experiments due 
to the sample roughness or the experimental treatments. 
Therefore, the quantum interference effect is different. 
Instead to analyzing this variation, we here confine our- 
selves to see the magnetic field dependence of Aa(B) and 
therefore depict the results for the normalized magneto- 
conductivity. The solid, dashed, and dashed-dot lines in 
Fig. 8 are the present calculations for the parameters 
(T, 5, L) = (0.12 K, 8 x 10~ 4 , (0.12 K, 1.5 x 10~ 3 , 

1/Ltm), and (7 K, 2 x 10~ 4 , 2/xm), respectively. These 
sets of parameters correspond/close to the conditions for 



the two experiments: the red filled circles and green 
filled squares from Ref. [l(| the up and down triangles 
and the diamonds from Ref. [T3- The impurity density 
rii = 1.15 x 10~ 3 a" 2 is the same as in our previous works 
which reproduce the zero-field electric conductivity^ and 
the Hall coefficient^ of the experimental results.— The in- 
tervalley scattering parameter is chosen as a = 0.7. Be- 
sides this, there is no additional adjustable free parame- 
ter. Clearly, the present calculation is in good agreement 
with the experimental measurements. 

Using the same diagrams as in Fig. 4, McCann et 
al. have also studied the magnetoconductivity adopt- 
ing constant parameters for disorders^ There are three 
scattering rates corresponding to three kinds of disorders 
in the formula for the magnetoresistivity. Here, because 
the charged impurities are considered as the predomi- 
nant scatters, except for the impurity density (the same 
as before for fitting the electric conductivity and the Hall 
coefficient of the experimental measurements^) , we have 
only the intcrvalley scattering as the parameter. 

The magnetoconductivity is due to the WL effect. By 
writing Eq. d3TJ) as SPi = ^; n ^-P/ n (similarly for 5P 2 ), 
we see that the contribution 8P n i(B) to the conductiv- 
ity from each Landau level of the center of mass of the 
Cooperon is, 



SPnl{B) cx 



B 



1 - X l (0) + 2di{2n + l)B' 



(32) 



For A°(0) = 1, 5P n o is independent of B. But the number 
of Landau levels in the summation depends on B. Since 
the g-integral of Cooperon is replaced by the summation, 
the singlet pseudospin channel gives rise to a positive 
contribution to Aa(B). For the pseudospin triplet, 1 — 
A'(0) > 0, and at small B « [1 - X l (0)]/2d l (2n + 1), 
5P n i(B) increases linearly with B. But the contribution 
to Sa(B) does not vanish at B = because the Landau 
levels become continuum and the g-integral is restored 
giving rise to a constant. The result of the summation 
at B = cancel with the corresponding term in c(0). 
Again, because of the summation over the discrete levels, 
the contribution to Aa(B) from the triplets is negative. 
Since the contributions from the pseudospin singlet and 
triplets vary differently as B varies, the final result for 
Aa(B) depends delicately on B. 

In the derivation given above, the single particle states 
were treated as the plane waves without taking into ac- 
count the Landau quantization. We here give the estima- 
tion of its validity. In the presence of the magnetic field 
B, the Landau level of the single Dirac-fermion particle 
is given b y 31 ' 32 
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FIG. 8: (color online) Normalized magnetoconductivity as 
function of the magnetic field B. The lines are the present 
calculations: solid, dashed, and dashed-dot lines are for the 
parameter sets (T,8,L) = (0.12 K, 8 x 10" 4 , IfJxa), (0.12 K, 
1.5 x 10~ 3 , l^rni), and (7 K, 2 x 10~ 4 , 2^m), respectively. 
The symbols are the experimental results: red filled circles 
(5 = 8 x 10~ 4 corresponding to n = 3.15 x 10 12 cm -2 ) and 
green filled squares [8 = 1.5 x 10~ 3 (n = 6.1 x 10 12 cm" 2 )], all 
at T — 0.12 K and L ~ l/j,m, are from Ref. I la : up triangles [T 
= 7K,8 = 2xl0~ 4 (n = 8x10" cm" 2 )], down triangles (T = 



E n = sgn(n)e 



(33) 



4 K, 8 = 2 x 10" 4 ), and diamonds [T = 0.26 K, S = 2.5 x 10" 4 
(n = 10 12 cm -2 )] all with L « 2^m are from Ref. Q]]. The 
magnetic field is given in unit of 10" 3 Tesla. 



doping 8 = 2 x 10" 4 and highest magnetic field B = 0.1 
Tesla considered here, the ratio reaches about 0.1. There- 
fore, the present calculation shown in Fig. 8 is valid. 



IV. SUMMARY 

In summary, on the basis of self-consistent Born ap- 
proximation, we have solved the Bethe-Salpeter ma- 
trix equations for Cooperon propagator of the Dirac 
fermions in graphene under the charged impurity scat- 
terings and weak external magnetic field. There are 
three non-vanishing components in the wavefunction of 
the Cooperon propagtor under the finite-range impurity 
scatterings. This feature is different from the zero-range 
one. The magnetoconductivity comes from the WL due 
to the quantum interference effect. The calculated mag- 
netoconductivity are in good agreement with the exper- 
imental measurements. 



with eo = vy // 2eB /he. The quantity eo is the over- 
all magnitude of the low level spacing. The above cal- 
culation of the magnetoconductivity is meaningful only 
when the ratio between eo and the Fermi energy Ep = 
(47T(5) 1 / 2 v/3 1 / 4 a is much less than unity. For the lowest 
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FIG. 9: (color online) The structure of a honeycomb lattice. 
There are two sites, a (black) and b (green), in each unit cell 
enclosed by the red lines. / is a vector from a to b. 
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APPENDIX 

In this appendix, we discuss the impurity scattering 
potential. In our previous work^ we have illustrated 
how the intravalley and intervalley scatterings vq(q) and 
i'i in Eq. |(3J) within the SCBA are determined from the 
microscopic electron- impurity interactions. In our nu- 
merical calculations, Vo(q) and v\ are set as respectively 
the values of leading terms in their expansions for the 
charged impurities. In the derivation, the difference of 
the a and b sites in the same unit cell of graphene lattice 
was neglected. As long as the long-wave length scatter- 
ings are considered, such a difference is negligible. It is 
true for the intravalley scatterings where the long-wave 
length scatterings are predominant momentum transfers 
of the Dirac fermions. While for the intervalley scatter- 
ings, the momentum transfers are finite and more careful 
treatment is needed. Here, taking into account the sub- 
lattice difference, we show that only the leading terms of 
vo(q) and v± need to be included in the effective poten- 
tials. 

We start with the Hamiltonian for electron-impurity 
interactions in graphene, 



rj)v ai (\fj - R\)rii(R) 



= y^2n a (q)v ai (q)m(-q), 



(34) 



qn 



where n a (fj) is the density operator of electrons at a 
(= a or 6) site of jth unit cell of the honeycomb lattice 
(Fig. 9), n.i(R) is the real space density distribution of 



impurities, and v a i(\rj — R\) is the impurity scattering 
potential. The Fourier component of the electron density 
n a(q) = J2k c l-qa c ku can be written as 

n »(q) ~ '(4-gal C feal + A-qa.1 C *oO) for q ~ 0, 
k 

!E '4- 9 'Q2 c fc"i for i ~ Q + <?' ~ Q> 
J2 c k- o'Qi c fc"2 for? Q + q Q, 
k 

where the subscripts 1 and 2 are the valley indices [see 
Fig. 1(a)] and k summation ^2 runs over a valley since 
we here consider only the low energy excitations. In the 
main text, the k summation ^ means £^ , and we here- 
after use the simple notation In Eq. (|34|) . the q 
summation runs over the infinitive momentum space. By 
folding the whole space into the first Brillouin zone, Hi 
can be written as 



kq 



(35) 



where k, q run over a valley in the Brillouin zone, and 



Vi(q) 



V d (q) 
Vj(q + Q) 



V (q-Q) 
Vd(q) 



with 

V d (q) 
V (q±Q) 

Mo) 

<f>b(Q) 



Mq) o 
o Mq) 



v d (q) = 



(36) 



Mq) o 
o Mq) 



o 

<t>°b(q±Q) 



Mq±Q) 
o 



J2 n i(-Q-Qn)V(q+Q n ) 

e^+Q^mi-q ~ Q n )V(q + Q n ). 



Here V(q) is the Fourier component of V a i(r), and Q n 
is the reciprocal lattice vector. There is a phase factor 
e i(q+Q n )-i m expansion of </>b(g) because of the posi- 
tion difference I between a and b sites. The expressions 
for 4>ab(o) are similar to 4> a ,b(q) except avoiding triple 
counting since <jf a b (q ± Q) imply the intervalley scatter- 
ings. [There are other two equivalent valleys for each 
valley indicated in Fig. QJa).] For example, the leading 
order 'm<%(q - Q) « <f>° a (-Q) is rn(Q)V(Q) where Q 
(with |Q| = 47r/3a, a ~ 2.4 A as the lattice constant) 
is a momentum difference between the nearest-neighbor 
Dirac points in the Brillouin zone. For this Q, there are 
other two vectors of the same magnitude differ from Q 
by two reciprocal lattice vectors [see Fig. Ufa)], respec- 
tively. These two terms should be excluded from the 
summation. Since the case of q << Q is under con- 
sideration, we approximate the off-diagonal potential as 
V (q±Q)^V (±Q). 

We have argueo~ that if the effective impurity scatter- 
ings are given as Eq. @, then the vertex T x (k,u)i,u)2) 
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can be expanded as Eq. ([5]). [The off-diagonal parts 
in Eq. @ are different from that in Ref. [l3| where 
cr was used instead of <j\ because the basis was = 
( c Li' c fc6i' c L2' c L 2 ) with reflected y-axis in valley 2. 
But the result for T x {k,ui\,u)2) is unchanged.] Here, we 



want to see how vo(q) and v\ in Eq. ((2]) are expected 
from the microscopic potentials given by Eq. (j5o) . The 
x-direction current vertex vY x {k, uji, L02) under the SCBA 
satisfies the following 4x4 matrix integral equation, 



I 

r x (fc, Wl ,o; 2 )=T3a x + ^^<^(fc-fcOG(^ (37) 

k' 

where (• • • ) means the average over the impurity distributions [Fig. He)]. Notice that the product 
G(k, iOi)T x (k, wi, w 2 )G(fc, C02) can be expanded in A x (k), 

G(k,uti)r x (k,LUi,LJ2)G(k,u>2) = y^^A*(k)Ljj>(k,uji,LJ2)yj'(k,uJi,L02), 

33' 

where the functions Ljji{k,u>x,u)2) have been defined in Ref. [HI. We then need to calculate the expectations (Vi(k — 
k^A^k'W^k' ~ k)) in Eq. (J37J). Firstly, we calculate (Vi(q)A%Vi(-q)) = A x [(<p a {q)<p b {-q)) - (MQ)atcM~Q)atc)- 
Using the expressions for (f> a ,b(Q)i we obtain 

(MQ)M-Q)) = Vm ]T V 2 (q + Q„)e- i W+^)- r « Vn % V 2 {q), 

n 

where the use of the fact that J2n^o \V(<1 + Qn)\ 2 e~ l ® n ' 1 = for q ~ has been made. Similarly, we get 
(MQ)atcM-QUc) = Vn t J2'v 2 (Q + Qn)e- l ^ + ^- r w VmV 2 {Q), 

n 

where J2 n ' means avoiding triple counting, and Q ■ I = because for I there is a Q orthogonal to it [see Figs. 1(a) 
and 9]. For A x (k') with j 7^ 0, we note 

A x 12 (k') = Al a (k){cos6±ia z 8m6), A%{k') = Ag(fc)(cos20 + ia z sin29), (38) 

where 9 is the angle between k and k . Since V(k' — k) depends on 9 through cos9, the integrals of the integrands 
with factor sin# or sin26 l vanish. (Vi(k — k')A x (k')Vi(k' — k)) for j ^ are then calculated as 

{Vi{k-k')A x 12 {k')Vi{k' - k)) -> A x 12 Vni ^V 2 {\k - k' + Q n \) cos9 « A x 12 V ni V 2 {\k! - fc|)cos<9 
(Vi{k- k')A%(k')Vi(k' - fc)) -> A%VmJ2V 2 (\k- k' + Q n |)e- 4( « +<3 " ) - r cos26» = AfFn;F 2 (|£' - fc|)cos26>. 

n 

By defining v 2 (q) = U 2 (g) and i; 2 = V 2 (Q), we obtain exactly the same equation as Eq. (11) in Ref. [l3|for determining 
yj(k,(jj±,W2). Thus, we have proved that only the leading terms of Vo(q) and v\ are necessarily taken into account in 
the current vertex corrections. 

On the other hand, under the SCBA, the self-energy is given by 

- ^]T(^(fc-fcOG(fc>Mfc'-fc)) 

k' 

« y X){E ^ - + Q ™) + E - Q")].9o(fc', w) + V 2 (k - k')g c (k', u)k ■ k'k ■ ut z } 

k' n n 

« y X^^* - ^) + ^ 2 (Q)]ff0(^', W) + U 2 (fc - *')Sc(fc', • fc'fc • #7i}, 

fc' 

I 

which is consistent with Eqs. © and Therefore, the theory is simplified with the effective potential given by 
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Eq. ([2]) with vo(q) and v\ denned above. 
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